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We present a computational study of the layered copper oxychloride (CuCl)LaNb2C>7 that has 
been recently proposed as a spin- 1/2 frustrated square lattice compound. Our results evidence an 
orbitally degenerate ground state for the reported tetragonal crystal structure and reveal a Jahn- 
Teller-type structural distortion. This distortion heavily changes the local environment of copper 
- CUO2CI2 plaquettes are formed instead of CUO2CI4 octahedra - and restores the single-orbital 
scenario typical for copper oxides and oxyhalides. The calculated distortion is consistent with the 
available diffraction data and the experimental results on the electric field gradients for the Cu and 
CI sites. The band structure suggests a complex three-dimensional spin model with the interactions 
up to the fourth neighbors. Despite the layered structure of (CuCl)LaNb207, the spin system 
has pronounced one-dimensional features. Yet, sizable interchain interactions lead to the strong 
frustration and likely cause the spin-gap behavior. Computational estimates of individual exchange 
couplings are in qualitative agreement with the experimental data. 

PACS numbers: 75.30.Et, 71.20.Ps, 71.70.Ej, 61.66.Fn 



I. INTRODUCTION 

Orbital ordering is one of the unusual and attractive 
phenomena in solid state physics. Specific occupation 
of orbitals plays a major role in electronic properties of 
numerous transition metal compounds. Orbital order- 
ing is generally accompanied by the Jahn- Teller effect, 
a structural distortion that lifts the orbital degeneracy.^ 
The Jahn- Teller distortion mainly affects the local envi- 
ronment of the transition metal cation, while the crys- 
tallographic unit cell is only slightly changed or even un- 
changed at all. The latter issue leads to certain difficul- 
ties in the structure analysis, since the distortion causes 
minor alterations of the diffraction patterns. Neverthe- 
less, modern experimental techniques of high-resolution 
neutron and synchrotron x-ray diffraction are usually 
able to resolve weak structural changes associated with 
the orbital ordering (see, e.g., Refs. I2I3I4I5[1 . 

Divalent copper is probably the most known transi- 
tion metal cation subjected to the Jahn- Teller effect. In 
a regular octahedral environment, the electronic configu- 
ration d 9 of Cu +2 leads to the orbital degeneracy that is 
usually lifted by an extremely strong tetragonal distor- 
tion. This distortion reduces the coordination number of 
copper and yields CUO4 plaquettes typical for the struc- 
tures of Cu +2 -containing oxides P High-symmetry struc- 
tures (e.g., perovskite structure) may constrain the Jahn- 
Teller distortion from the plaquettes formation. Never- 
theless, a notable tetragonal distortion takes place, the 
orbital ordering is established, and peculiar electronic 
properties emerge. Thus, perovskite- type copper fluo- 
rides KCUF3 and K2CUF4 are known as first examples 
of the cooperative Jahn- Teller distortion, giving rise to 
ferromagnetic ( FM) interactions induced by the specific 
orbital ordering! 1 ! 6 ! 7 ! 

Despite the large magnitude in the Jahn- Teller distor- 



tion in copper compounds, the identification of the dis- 
tortion pattern may be quite problematic. This is the 
case for layered copper oxyhalides (CuX)LaM 2 07 with 
X = CI, Br and M = Nb, Ta. These compounds have 
composite structures built by [LaM^Or] perovskit e-ty pe 
blocks and [CuX] rocksalt-type layers (see Fig. JTJPS In 
the following, we will focus on one of these compounds, 
(CuCl)LaNb 2 07, that recently drew atte ntion d ue to its 
unusual and puzzling magnetic properties! 10 ! 11 ! 

Initially, the crystal structure of (CuCl)LaNb207 was 
refined in the tetragonal space group PA/mmm with 
the CI atom located in the special 16 position (0,0, |)P 
In this structure (further referred as regular), copper 
has a squeezed octahedral coordination with two short 
Cu-0 bonds [d(Cu-O) = 1.97 A] and four long Cu-Cl 
bonds [<f (Cu-Cl) = 2.74 A], see Fig. [I] This type of 
local environment is quite unusual for copper oxychlo- 
rides: normally, these compounds reveal the Jahn- Teller 
effect and the square-planar CuCU or CUO2CI2 coordina- 
tion with d(Cu-Cl)= 2.3 - 2.4 AP Moreover, the refined 
Debye-Waller factor for the CI atom is extremely high 
{U\ so — 0.13 A 2 ) suggesting a shift of CI away from the 
lb positionP Caruntu et aZ™ proposed a new structural 
model with the CI atoms randomly occupying one quarter 
of 4m sites (x,0, |) with x = 0.136. This model yields 
two short (about 2.4 A) and two long (about 3.15 A) 
Cu-Cl distances consistent with the crystal chemistry 
of copper oxychloridesP Caruntu et aZP^l tentatively as- 
cribed the distortion to the Jahn- Teller effect of Cu +2 . 
However, they failed to observe any superstructure reflec- 
tions that could arise due to the cooperative Jahn- Teller 
distortion, similar to KCuF 3 and K 2 CuF 4 P^I 

The structure of (CuCl)LaNb20 7 was further studied 
by means of nuclear magnetic resonance (NMR) and nu- 
clear quadrupole resonance (NQR) measurements.^ The 
spectra evidenced the lack of the tetragonal symmetry 
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FIG. 1: (Color online) Regular (tetragonal) crystal structure 
of (CuCl)LaNb2 07: stacking of perovskite-type [LaNb2 07] 
blocks and rocksalt-type [CuCl] layers (left panel), the single 
[CuCl] layer (upper right panel), and the CUO2CI4 squeezed 
octahedron (bottom right panel). The regular structure of 
the [CuCl] layer gives rise to the frustrated square lattice spin 
model with the competing nearest-neighbor and next-nearest- 
neighbor couplings Ji r and .hr, respectively. 



and revealed single sites of Cu, CI, and La atoms, hence 
suggesting the ordering of the CI atoms, at least on the 
local scale. Additionally, Yoshida et al. 13 performed elec- 
tron diffraction studies and found superstructure reflec- 
tions that unambiguously confirmed the ordering of the 
CI atoms and the resulting cooperative distortion of the 
copper polyhedra. However, the specific ordering pattern 
remains unclear. To refine the superstructure, one has to 
use the superstructure reflection s that a re absent in x-ray 
and neutron diffraction patterns! 8 * 12 * 14 * The respective re- 
flections are revealed by electron diffraction, but their in- 
tensities are strongly affected by multiple scattering and 
hence unsuitable for the refinement.^ 

The magnetic properties of (CuCl)LaNb207 are un- 
usual and lack a clear microscopic interpretation. The 
studies reveal a spin gap behavior and a singlet ground 
state that are incompatible with the frustrated square 
lattice (Ji r — J2r) m odel, as p roposed by the regular crys- 
tal structure in Fig. [lp^ 11 * 16 * Inelastic neutron scattering 
(INS) data are even more puzzling and suggest long-range 
interactions (between Cu atoms separated for about 9 A) 
to be relevant for (CuCl)LaNb207p^ Such an unusual 
magnetic behavior could be caused by a non-trivial or- 
bital state of copper and intricate superexchange path- 
ways, emerging due to a specific cooperative distortion. 
Thus, unraveling this distortion is a key to understand 
the magnetic properties of (CuCl)LaNb 2 07. 

Presently, there are two reports that consider pos- 
sible ordering patterns of the CI atoms and attempt 
to relate these patterns to the exchange interactions in 
(CuCl)LaNb207. Whangbo and Daff^ used extended 



Hiickel calculations to study qualitatively exchange cou- 
plings in several ordered structures. They proposed a 
model of ring clusters with a number of incquivalent Cu 
and CI sites that is in contradiction to the NMR and 
NQR results.^ Yoshida et al. 13 employed an even more 
empirical approach and searched for an ordering pattern 
with "explicit" copper dimers (i.e., two Cu atoms con- 
nected by the double bridge of short Cu-Cl bonds) . Such 
dimers can probably account for the spin-gap behavior, 
but they are unable to explain the long-range interac- 
tions, as proposed by the INS data. Neither of the two 
studies considered the underlying reasons for the struc- 
tural distortion or proposed numerical estimates of the 
exchange couplings. 

The physics of (CuCl)LaNb 2 07 is quite complex and 
includes both structural (cooperative distortion) and 
magnetic (exchange interactions) aspects. Unfortunately, 
experimental techniques up to now fail to determine the 
superstructure, hence complicating the analysis of the ex- 
change couplings. The latter is non-trivial itself due to 
the presence of long-range interactions. In this work, we 
attack the problem by using full-potential band struc- 
ture calculations based on the density functional the- 
ory (DFT), since computational methods are known as 
an effective tool for solving complex problems in solid 
state physics. The main advantage of the computa- 
tional approach is the possibility to estimate numeri- 
cally the relevant parameters such as the total energy 
and the exchange couplings. The evaluation of the to- 
tal energy is crucial for resolving the structural prob- 
lems in case conventional diffraction techniques experi- 
ence difficulties in a precise structure determination. For 
example, the local density approximation (LDA) level 
of DFT was successfully applied for studying the struc- 
tures of lithium monoboride, 18 while proper treatment 
of correlation effects allowed to derive Jahn- Teller dis- 
tortion and to explain the resulting FM interactions in 
transition-metal fluorides KCUF3 (Refs. [19] and l20jl and 
Cs2AgF4 (Ref. 21). Quite recently, band structure cal- 
culations predicted ferroelect ric dis tortions in unconven- 
tional manganite multiferroics^^ followed by the exper- 
imental confirmation!^ 

The outline of the paper is as follows. We start with 
the discussion of methodological aspects in Sec. ITT] In 



Sec. Ill we present the computational results for the reg- 
ular structure of (CuCl)LaNb 2 07. We find orbital degen- 
eracy and show that the regular structure is in qualita- 
tive contradiction to the experimental magnetic behavior. 
In Sec. |IV[ we derive the relevant structural distortion, 
evaluate exchange integrals, and propose a realistic spin 
model of (CuCl)LaNb 2 7 . In Sec. |V) we compare our 
results to the experimental findings, discuss the resulting 
spin model, and demonstrate that our scenario provides 
consistent interpretation of the experimental data. We 
also give an outlook for further experimental studies and 
conclude the section with a short summary. 



3 



II. METHODS 

Scalar-relativistic band structure calculations were 
performed using a full-potential local-orbital scheme 
(FPLO, versions 7.00-27 and 8.50-32)P and the 
exchange-correlation potential by Perdew and Wang™ 
LDA calculations were done for the smallest possible unit 
cells employing the actual symmetry of the atomic posi- 
tions, i.e., a 12-atom tetragonal unit cell for the regular 
structure (a x a x c, space group PA/mmm) and a 24- 
atom orthorhombic unit cell for the distorted structure 
(a x 2a x c, space group Pbmm). For the local spin 
density approximation (LSDA)+C/ calculations with dif- 
ferent types of spin ordering, a number of monoclinic 48- 
and 72-atom supercells were used (space groups Pm or 
P2/m with c being the monoclinic axis, see Sec. IV and 
Fig. [8] for details). The structure optimization was car- 
ried out in the 48-atom 2a x 2a x c supercell, as suggested 
by the electron diffraction dataP^The triclinic (PI) sym- 
metry of this supercell allowed for free relaxation of all 48 
atoms. Different unit cells employed different k meshes 
with at least 72, 192, 730, and 2925 k points in the ir- 
reducible part of the first Brillouin zone for the 72-, 48-, 
24-, and 12-atom cells, respectively. The convergence 
with respect to the k mesh was carefully checked. 

Two sets of structural information are available for 
(CuCl)LaNb207. In Ref. the structure was refined 
from the x-ray data, while in Ref. [12] the neutron data 
were used. The latter results are claimed to be more ac- 
curate, because neutrons are more sensitive to the posi- 
tions of light atomsP2l Additionally, the interatomic dis- 
tances yield more reasonable formal valences, as calcu- 
lated with the empirical bond valence sum rulesP^ Yet, 
the authors of Ref. |T2] note that the resulting Cu-0 dis- 
tance (1.84 A) is unusually short as compared to other 
copper oxides and oxychloridesP Keeping in mind this 
ambiguity, we performed calculations for both sets of 
the atomic coordinates. The results for the stabilization 
energies and the exchange couplings are similar within 
5 — 10%, which is likely comparable to the actual accu- 
racy of our calculations. This outcome is quite natural, 
because the difference between the two structures deals 
with the atomic positions along the c axis (the positions 
in the ab plane are fixed by the symmetry), whereas the 
main physics takes place in the ab plane. Note, how- 
ever, that the structure refined from the neutron data 
has lower energy, though both structures are not fully re- 
laxed and reveal sizable forces (0.3 — 0.4 eV/a.u.) along 
the c axis for Nb and O atoms. The forces on the respec- 
tive atoms have opposite signs, hence the relaxed atomic 
positions lie between those suggested by x-rays and neu- 
trons. 

The general computational procedure for both regu- 
lar and distorted structures is as follows. We start with 
LDA calculations, focus on the states close to the Fermi 
level and construct an effective model for these states. 
The model yields the estimates for all the antiferromag- 
netic (AFM) couplings in the systems under investiga- 



tion, hence facilitating the choice of relevant interac- 
tions. We also calculate Wannier functions in order to 
get a direct, graphical representation of the relevant or- 
bitals. The Wannier functions are computed in a max- 
imally localized fashion 2 '' using the internal procedure 
of FPLO8.50. Then we proceed to LSDA+C/ calcula- 
tions and get independent estimates of the leading ex- 
change couplings. Additionally, we use LSDA+C/ for the 
structure optimization, because LDA is known to conceal 
Jahn-Teller distortions in transition metal compo unds 
due to improper treatment of the correlation effectsP^D 

In the LDA band structure, the states close to the 
Fermi level are relevant for the possible orbital degen- 
eracy and for the magnetic interactions. The respec- 
tive bands are analyzed within a tight-binding (TB) 
model, and the resulting hoppings (t) are introduced 
to an extended Hubbard model with the effective on- 
site Coulomb repulsion U c g. For both regular and dis- 
torted structures, the strongly correlated limit t <C J7 e ff 
is realized, hence the Hubbard model can be reduced to 
the Heisenberg model for the low-lying excitations un- 
der the assumption of the half-filling. The AFM con- 
tributions to the exchange integrals are estimated as 
jAfm _ 4^/i/ eSj an d a n the possible superexchangc 
pathways are analyzed. The latter feature is especially 
important, since long-range interactions are expected in 
(CuCl)LaNb 2 7 PThe precise value of the U e g param- 
eter for copper oxyhalides is unknown, although a re- 
cent study of CuCl 2 and CuCl2 • 2H 2 suggests U c s = 
4 eV) 28 * 29 ! In the present work, we employ the represen- 
tative value U e g = 4.5 eV, as commonly used for copper 
oxides (see, e.g., Refs. l30l31l32l33p . Note that the values 
of J AFM merely scale with U e g, and the ratios of J AFM , 
which arc the relevant quantity for the ground state, are 
not affected by the ambiguity of the choice of U e g. 

The LDA results can be effectively used to construct 
models on top of DFT. However, it is also helpful to ac- 
count for the correlation effects in the course of the self- 
consistent calculations. This opportunity is provided by 
the LSDA+C/ method that treats the correlation effects 
in the mean-field approximation. LSDA+C/ is known to 
yield reasonable estimates of different quantities [energy 
gaps, electric field gradients (EFGs), and exchange in- 
tegrals] in strongly correlated electronic systems and, in 
particular, to tr eat correctly the cooperative Jahn- Teller 
distortions! 19 * 21 ! The problematic point of the LSDA+C/ 
approach deals with the choice of the Ud and Jd val- 
ues that account for the on-site Coulomb repulsion and 
exchange, respectively. The optimal values of these pa- 
rameters are known to depend on the particular com- 
putational method and on the objective quantities as 
well as on the specific chemical and st ructura l features 
of the compound under investigationPSESEI] Unfortu- 
nately, there are no established LSDA+C/ parameters for 
copper oxy chlorides, although the recent application of 
U d = 6.0 - 8.5 eV to C uCl 2 and CuCl 2 • 2H 2 can be 
used as a reference ! 28 * 29 * In our calculations, we employ 
several representative and physically reasonable values 



4 



of Ud over a wide range, namely, 3.5, 5.5, and 7.5 eV.^] 
Below, we will show that all these values lead to quali- 
tatively similar behavior, although the details of the mi- 
croscopic scenario depend on the Ud value used. The 
exchange parameter Jd is fixed at 1 eV, since it is known 
to have minor influence on the results. 

One should keep in mind that the LSDA+C/ Coulomb 
repulsion parameter Ud is usually different from the ef- 
fective on-site Coulomb repulsion U c g, as employed in the 
model analysis above. The latter potential corresponds 
to the antibonding pda* LDA bands with the contribu- 
tions from Cu, CI, and O orbitals (see Fig. [2]), while Ud 
is applied to the Cu 3d orbitals only. In the framework 
of the one-band Hubbard model, the antiferromagnetic 
superexchange couplings should scale as 1/U c g, but this 
scaling is not necessarily applicable to the LSDA+[7 re- 
sults for the total exchange and to their dependence on 
U d . 



III. REGULAR STRUCTURE 

The LDA density of states (DOS) for the regular 
structure of (CuCl)LaNb 2 07 is shown in Fig. [2j The 
energy spec trum is similar to that in other copper 
oxy chlorides ESES1 The low- lying valence bands are pre- 
dominantly formed by oxygen orbitals, while at higher 
energies copper and chlorine states become prominent. 
The states near the Fermi are composed of copper or- 
bitals, with sizable contributions from oxygen and chlo- 
rine. Conduction bands are formed by niobium and oxy- 
gen orbitals. The spectrum is consistent with the intu- 
itive ionic picture, suggesting the oxidation states of +2 
and +5 for Cu (3d 9 ) and Nb (4<f°), respectively. Thus, 
peculiar magnetic properties originate from the [CuCl] 
layers, while the Nb-containing [LaNb 2 Oy] blocks are 
non-magnetic and insulating. 

The energy spectrum is metallic in contradiction to 
the observed green color of (CuCl)LaNb207P The unre- 
alistic metallicity is a typical failure of LDA due to the 
underestimate of the correlation effects. The application 
of LSDA+J7 readily restores an energy gap of 1.0 — 1.2 
eV for both the regular and distorted structures. Such 
a gap is still insufficient to account for the green color 
of the compound. However, one should be aware that 
the straightforward mapping of either the LDA or the 
LSDA+J7 single-particle spectrum onto the optical prop- 
erties may be misleading. Many-body effects are some- 
times crucial for proper evaluation of the excitation en- 
ergies, as recently shown for lithium niobate, LiNb03. 37 
We believe that the optical properties of (CuCl)LaNb207 
require a careful study on the many-body level of theory, 
and such a study lies beyond the scope of the present 
work. In the following, we restrict ourselves to the discus- 
sion of ground-state properties (optimized crystal struc- 
ture) and low- lying excitations (magnetic interactions). 
Both problems deal solely with the valence bands, while 
the conduction bands lie 1 — 2 eV higher in energy hence 
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FIG. 2: (Color online) LDA density of states (DOS) for the 
regular (tetragonal) structure of (CuCl)LaNb207. The Fermi 
level is at zero energy, the states close to the Fermi level are 
formed by Cu, CI, and O orbitals. The inset shows a part of 
the orbital resolved DOS for copper with two orbitals (d 3z 2_ r 2 
and d x 2_ y 2) contributing to the upper valence bands. 



bearing little influence on the objectives of our work. 

The orbital resolved DOS shows that the bands near 
the Fermi level are mainly formed by Cu d 3z 2_ r 2 orbital 
with a contribution from the d x 2_ y 2 orbital (see the inset 
of Fig. [2| . The same picture emerges from the band struc- 
ture plot. There are two distinct copper bands above 
— 1 eV, and these bands are formed by the d^ z 2_ r 2 and 
d x 2_ y 2 orbitals (Fig. |3|. Crystal field considerations sug- 
gest that the orbital with the highest energy should point 
to the nearest ligands. Indeed, the d 3z 2_ r 2 orbital points 
to oxygen atoms with short Cu-0 separations of about 
1.85 A. The in- plane Cu-Cl separations are much longer 
(2.74 A), but the chlorine orbitals are spatially more ex- 
tended as compared to the oxygen orbitals. Thus, both 
e g orbitals of copper have similar energies, and the LDA 
ground state is orbitally degenerate. Correlations induce 
an orbital ordering within the regular structure, but - as 
we will show below such orbital ordering is unable to 
account for the magnetic properties of (CuCl)LaNb20y. 

The d^ z 2_ r 2 band can be analyzed within a simple 
TB model, including nearest-neighbor (t\ r ) and next- 
nearest-neighbor (t2r) hoppings (the subscript r refers 

<i 2 _ 2 

to the regular structure). We find t^, r = —3 meV 

and t~2r r — 183 meV. The fitting of the d x 2_ y 2 band 
is more complicated due to its strong hybridization with 
the chlorine orbitals. However, one can roughly estimate 

t\ r v = —44 meV and t^ r ~ v — 172 meV from the 
r — X — M part of the band structure. Thus, both or- 
bitals yield the t 2r ^> t\ r scenario and suggest the strong 
AFM interaction J 2 A r FM of about 300 - 350 K. 

Electron correlations are able to lift the orbital de- 
generacy without involving a structural distortion, as 
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FIG. 3: (Color online) LDA bands for the regular (tetrag- 
onal) crystal structure of (CuCl)LaNb2 07. Lines show the 
bands, while dots indicate the contributions of the d 3z 2_ r 2 
and d x 2_ y 2 orbitals. The Fermi level is at zero energy. The 
notation of k points is as follows: r(0,0,0), A (0.5,0, 0), 
A/(0.5, 0.5,0), Z(0, 0,0.5), R(0.5, 0.5, 0.5), and A(0.5, 0, 0.5), 
where the coordinates are given along the k x , k y , and k z axes 
in units of the respective reciprocal lattice parameters 2n/a, 
2n/a, and 2n/c. 

recently shown for CuSb2C>6.^ In this compound, the 
LDA band structure reveals similar contributions of the 
d-i z 2_ r 2 and d x 2_ y 2 orbitals to the states near the Fermi 
level. The introduction of electron correlations within 
LSDA+J7 enables to stabilize the half- filling for either 
of the orbitals, but the d^ z 2_ r 2 -type orbital ordering 
has lower energy and yields the realistic physical pic- 
ture. The case of (CuCl)LaNb 2 7 is different. LSDA+C/ 
leads to the half-filling of the d 3z 2_ r 2 orbitals (the state 
with the d x 2_ y 2-type ordering can not be stabilized), 
but the resulting scenario contradicts the experimen- 
tal results. Applying Ud — 3.5 — 7.5 eV, we find the 
J\ r values in the range from —50 K to —10 K and 
J2r = 240 — 360 K. These numbers are in reasonable 
agreement with the TB estimates but in conflict with 
the experimental energy scale. Experimental values of 
the Curie- Weiss temperature (9 = 9.6 K)G^and the satu- 
ration field (fioH s — 30.1 T}^are quite small. Assuming 
the frustrated square lattice model with J lr = —30 K 
and J 2r = 300 K, we find 9 = J lr + J 2r = 270 rPH and 
fx H s = (2Ji r + 4J 2 r)W(.9Ms) ~ 840 T, 39 in obvious 
contradiction to the experimental results. Thus, correla- 
tions alone are unable to account for the orbital ordering 
and the resulting spin physics of (CuCl)LaNb2C>7. Below, 
we demonstrate that one has to consider lattice degrees 
of freedom as well. 



IV. REALISTIC STRUCTURE 

A. Structural distortion 

To derive the structural distortion in (CuCl)LaNb2C>7, 
we use the supercell with doubled a and b parameters, as 
suggested by the electron diffraction studies P^l We focus 



on the forces in the ab plane and neglect the forces along 
the c axis. In Sec. [TTJ, we have noticed the sizable forces 
along the c axis for Nb and O atoms due to the uncer- 
tainty of the experimental structure determination. Such 
forces are present in both the regular and the distorted 
structures. Our calculations are performed for two dif- 
ferent sets of structural data showing forces of opposite 
signs. The two sets are opposite with respect to the fully 
relaxed one and yield essentially similar results. Hence, 
the precise positions of the atoms along the c axis bear 
little influence on the structure in the ab plane and on 
the resulting magnetic interactions. 

The regular structure does not reveal any sizable forces 
in the ab plane. However, the displacement of the CI 
atoms away from the high-symmetry position leads to 
the forces that tend to enhance the distortion. The su- 
percell includes four independent CI atoms, but the forces 
on these atoms are found to be similar. This indicates 
the cooperative character of the distortion. The shifts 
of the CI atoms give rise to two short and two long Cu- 
Cl bonds. The resulting distortion pattern is shown in 
Fig. Ptl The CI atoms occupy the (x, 0, |) sites in an or- 
derea manner to yield opposite directions for the short 
Cu-Cl bonds on each Cu site. The structural transfor- 
mation is quantified by a single parameter, the x value 
in (x,0, i) (in the following, we assume x to be the co- 
ordinate of the shifted CI atoms in the initial, tetragonal 
unit cell). The change in the total energy in the course 
of the distortion is visualized in Fig. [5] 

Both LDA and LSDA+J7 predict the stabilization of 
the distorted structure with respect to the regular one. 
The relaxed positions of the CI atoms are nearly the 
same in all the calculations, while the stabilization en- 
ergy (AE) shows a pronounced dependence on the com- 
putational method and on Ud- The AE value is the 
energy gain due to the formation of two short Cu-Cl 
bonds in the distorted structure rather than four longer 
bonds in the regular structure (see the upper part of 
Fig. [Sj. Therefore, AE depends on the hybridization 
between Cu and CI orbitals, the latter being controlled 
by the Ud value. The application of the Ud potential 
shifts filled copper states to lower energy, hence pro- 
viding stronger hybridization and increasing AE at low 
Ud (3.5 eV). Higher values of Ud further shift the cop- 
per states to lower energy, reduce the hybridization, and 
AE is decreased. Yet, even the smallest (LDA) value 
(AE = 0.15 eV/f.u.) is sufficient to stabilize the distor- 
tion with respect to thermal fluctuations (at room tem- 
perature, thermal energy amounts to about 0.025 eV). 
This finding contrasts to the studies of KCUF3 (Ref. 
and Cs2AgF4 (Ref. [3TJ, where the distortion is stabilized 
within LSDA+f/ only.' In (CuCl)LaNb 2 7 , both the lat- 
tice and the correlations favor the distorted structure. 

Our relaxed structure is different from one suggested 
by Yoshida et al. 13 Basically, an ordering similar to ours 
was considered by Whangbo and Daif^ though they did 
not optimize the positions of the CI atoms, did not find 
the opportunity for the spin-gap behavior and, finally, 



FIG. 4: (Color online) Structural distortion in the [CuCl] layers: the forces on the CI atoms (left panel), magnetic interactions 
in the distorted layer (middle panel), and the resulting spin model (right panel). Smaller and larger spheres show the Cu and 
CI atoms, respectively; <p measures the Cu-Cl-Cu angle. In the left panel, the square marks the supercell employed in the 
structure optimization, while in the middle panel the rectangle shows the crystallographic unit cell of the distorted structure. 
Thick lines in the middle panel are the short Cu-Cl bonds, indicating the CUO2CI2 plaquettes arranged perpendicular to the 
layer. In the right panel, only the leading interactions are shown. The solid, dotted, and dashed lines denote Ji, J2, and J4, 
respectively. 



0.0 



3 -0.1 

> 

^ -0.2 
-0.3 




LSDA+l/, U, 

—a— 3.5 eV 
— •— 5.5 eV 

— o— 7.5 eV 



^0-0-0-0-0- 







0.05 



0.10 

x(Cl) 



0.15 



0.20 



FIG. 5: (Color online) The transformation of the copper 
local environment (upper part) and the change in the to- 
tal energy (bottom part) along the structural distortion in 
(CuCl)LaNb207. The total energy of the regular (tetragonal) 
structure is set to zero. The upper axis shows the evolution 
of the shorter Cu-Cl distance. In the bottom axis, x indi- 
cates the coordinate of the (x,0, |) position of the CI atoms. 
The x value is measured in units of a (the lattice parame- 
ter of the regular structure). The dashed line shows the x 
coordinate of the experimental four-fold position of CI with 
random occupation, as determined by the refinement of the 
neutron dataP^" In the upper panel, the shaded region marks 
the CUO2CI2 plaquette. 



discarded the model. Neither of the two groups sup- 
ported the choice of the model by total energy calcula- 
tions and proved the stability of the distorted structure 
with respect to the regular one. In contrast to the pre- 
vious studies, our proposition for the structural model is 
well-justified and confirmed by the sizable stabilization 



energies on different levels of theory (Fig. [5j. Below, we 
will demonstrate that our model is also consistent with 
the available experimental data. 

First of all, the shift of the CI atoms to the (a;,0, \) 
position is in agreement with the neutron diffraction 
results Despite the lack of the superstructure reflec- 
tions in the neutron data, it was possible to refine the x 
value by assuming random occupation of four equivalent 
(x, 0, |) sites. The resulting x ~ 0.136 (dashed verti- 
cal line in Fig. |5| is in reasonable agreement with the 
positions of the energy minima (x = 0.16 — 0.18). The 
optimized structure shows single sites for the Cu, CI, and 
La atoms in agreement with the NMR and NQR dataP^l 
Yet, the a x 2a x c unit cell (see Fig. |4]) seems to contradict 
the electron diffraction results, suggesting the doubling 
of both a and b parameters.^ However, one should keep 
in mind that in the regular structure the a and b axes 
are indistinguishable, hence the CI atoms can be equally 
shifted to either (x,0, |) or (0,x, \). Then, one expects 
the 90° twinning with twins leading to two different sets 
of superstructure reflections in the electron diffraction 
patterns and to the apparent doubling along both the a 
and b axes. Thus, the twinning can explain the seeming 
controversy between the structural model and the elec- 
tron diffraction data, although a careful high-resolution 
electron microscopy (HREM) study is desirable to resolve 
this issue. 

We find a further experimental confirmation for the 
proposed structural model by considering EFG's. The 
EFG's for the Cu, CI, and La sites were measured exper- 
imentally by NMR on aligned powder samples.- 1 -"* There- 
fore, several characteristics can be determined: (i) the 
leading principal component of the tensor (V zz ); (ii) the 
asymmetry parameter [77 = (V xx — V vy )/V zz ); and (iii) the 
orientation of the z axis of the tensor with respect to the 
c axis of the crystal structure. The EFG's are a sensitive 
tool to detect the local symmetry of an atomic position. 
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TABLE I: Computational estimates of the electric field 
gradients for the distorted and regular structures of 
(CuCl)LaNb207: leading principal components of the ten- 
sors V zz (in 10 21 V/m 2 ), asymmetries r), and orientations of 
the tensor (2 vs. c). Ud is the Coulomb repulsion parameter 
of LSDA+17. The bottom part of the table lists the experi- 
mental results from Ref. 1131 
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"Note that the NMR data do not allow to determine the sign 
of Vzz- The values are listed according to Ref. 1131 though the 
assignment is basically arbitrary. 



In tetragonal structures, the positions on the four-fold 
axes should have V xx = V yy , hence r\ = 0. The non-zero 
experimental rj value for the Cu, CI, and La sites (see the 
last line of Table [I]) is a strong evidence for the structural 
distortion in (CuCl)LaNb 2 7 P 

The calculated EFG parameters for the regular and 
distorted structures are listed in Table [i] along with the 
experimental results. The regular structure suggests 
77 = for all the three sites and z || c for the Cu site 
in clear contradiction to the experiment. Yet, the dis- 
torted structure provides reasonable r\ values for the Cu 
and CI sites. For the Cu site, V zz and 77 strongly depend 
on Ud, since the on-site correlations modify the electronic 
distribution around this site. For the CI atoms, the de- 
pendence on Ud is weak, while the results for the La site 
are almost independent on Ud, because the La and Cu 
sites are far apart. At Ud = 5.5 eV and 7.5 eV, the EFG's 
on the Cu and CI sites are in good agreement with the 
experiment. The orientation of the EFG tensors is also 
reproduced. For the La site, the calculations yield the 
reasonable value of V zz and the correct orientation of the 
tensor. However, the calculated asymmetry parameter is 
well below the experimental one. 

To understand the discrepancy for the asymmetry of 
the La site EFG, one should recall that the low asym- 
metry implies the weak distortion of the local envi- 
ronment with respect to the regular structure (for the 
regular structure 77 — 0, see Table Our structural 
model suggests the distortion of the [CuCl] layers, while 
the [LaNb207] blocks retain the tetragonal symmetry. 
Therefore, the low calculated 77 value is reasonable. The 
experimental 77 = 0.70 implies the distortion of the 
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FIG. 6: (Color online) LDA bands for the distorted struc- 
ture of (CuCl)LaNb207 (left part) and the respective den- 
sity of states for copper atoms (right part). The Fermi level 
is at zero energy. The CI atoms occupy the (x,0, |) po- 
sition with x — 0.16. In the left panel, the thin (light) 
lines denote the LDA bands, while the thick (dark) lines 
show the fit of the tight-binding model. In the right panel, 
the solid line indicates the total density of states for copper 
atoms, and the solid filling marks the contribution of the in- 
plaquette orbitals. The notation of k points is as follows: 
r(0,0,0), X(0.5,0,0), M(0.5, 0.5,0), Y (0,0.5,0), Z(0, 0,0.5), 
4(0,0.5,0.5), #(0.5,0.5,0.5), and T(0.5, 0,0.5) (the coordi- 
nates are given along the k x ,k y , and k z axes in units of the 
respective reciprocal lattice parameters 2n/a, n/a, and 2n/c). 



[LaNb207] blocks. Such a distortion may take place due 
to a tilting of the NbOe octahedra. The t ilting distortions 
are typical for perovskite-type structures and one can 
expect this type of the distortion in (CuCl)LaNb20 7 P3 
However, the band structure calculations do not evidence 
any changes in the [LaNb20 7 ] blocks, since the ab com- 
ponents of the forces on the Nb and O atoms are negli- 
gible (below 0.01 eV/a.u.) The high-symmetry structure 
may correspond to a local energy maximum, where the 
forces should vanish. However, the shifts of the oxygen 
atoms away from the high-symmetry positions yield siz- 
able forces that tend to restore the high-symmetry struc- 
ture. Thus, there is no evidence for the distortion of the 
[LaNb207] blocks. Possibly, one has to use a larger su- 
percell in order to observe this distortion. However, the 
electron diffraction datcP do not suggest a larger super- 
cell leaving the distortion of the [LaNb207] blocks as an 
open question. 

Finally, one finds support for the distorted structure 
by considering the crystal chemistry of copper com- 
pounds. The distortion leads to the planar CUO2CI2 
environment of copper (see the upper part of Fig. [5]) 
with two Cu-Cl bonds of 2.3 - 2.4 A. Two other Cu-Cl 
distances are extended up to 3.2 A and become non- 
bonding. The CUO2CI2 planar coordination is typical 
for copper oxychlorides, and the CI atoms are normally 
located in the opposite corners of the O2CI2 rectangle 
similar to our model (Fig. [4|. 6 Therefore, the distortion 
in (CuCl)LaNb207 is also in accord with the empirical 
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structural expectations. 

Basically, the CUO2CI2 unit can be considered as 
an analog of the conventional CUO4 plaquette. The 
CUO2CI2 plaquette is geometrically distorted due to the 
different Cu-0 and Cu-Cl bond lengths. However, its 
electronic structure is similar to that of the CUO4 square. 
The half-filled orbital lies in the CUO2CI2 plane, and the 
lobes of this orbital point to both the O and CI atoms, 
as illustrated by the respective Wannier function (see the 
left panel of Fig. |9j. The Wannier function also shows 
sizable contributions from u-type O and CI orbitals. The 
position of the copper orbital is in agreement with intu- 
itive crystal-field considerations. LDA reveals two bands 
at the Fermi level, and these bands are predominantly 
formed by the in-plaquette orbitals of two copper atoms 
(see Fig. [6]). Thus, the structural distortion lifts the or- 
bital degeneracy and leads to the single-orbital ground 
state typical for most of the Cu +2 compounds. 



B. Exchange couplings 

In this section, we evaluate the leading ex- 
change couplings for the distorted crystal structure of 
(CuCl)LaNb207 and derive a realistic spin model for this 
compound. Below, we use two different (x, 0, |) positions 
of the CI atoms: x — 0.16 as a representative position 
for the optimized structure (see Fig. |5| and x — 0.136 
as the coordinate suggested by the neutron diffraction 
experiment We also consider general trends for the 
evolution of the exchange couplings along the shift of 
the CI atom. 

The distortion lifts the tetragonal symmetry and leads 
to a large number of inequivalent superexchange path- 
ways (Fig. [4|. The distorted layer can be considered as 
a system of corner-sharing CUO2CI2 plaquettes forming 
chains along the b axis. Then, the single nearest-neighbor 
interaction of the regular structure (Ji r , see Fig. [T]) is 
split into two interactions: J\ runs along the chains and 
J[ runs perpendicular to the chains. The next-nearest- 
neighbor coupling J2 is unique, similar to the regular 
structure. Third-neighbor interactions are again split 
into two interactions: J3 runs along the chains and J3 
runs perpendicular to the chains. Finally, there are three 
inequivalent fourth-neighbor couplings J4, J' A , and J'l . In 
the following, we select the relevant interactions by ana- 
lyzing the LDA band structure. Then, we also evaluate 
these interactions via the LSDA+[/ calculations. 

The LDA band structure is shown in Fig.[6j Two bands 
at the Fermi level correspond to two Cu atoms per dou- 
bled unit cell in contrast to Fig. [3] with two bands, origi- 
nating from two orbitals of the same Cu atom in the orig- 
inal unit cell. The bands in Fig. [6] are fitted with a TB 
model. The parameters of this model (hoppings) are cal- 
culated as overlap integrals of the Wannier functions cen- 
tered on copper sites. The resulting numbers are listed 
in Table [II] The largest hopping runs along the chains of 
the CUO2CI2 plaquettes and leads to the AFM interac- 




FIG. 7: The evolution of the leading hoppings parameters 
within the [CuCl] layer upon shifting the CI atom. The x 
value is measured in units of a (the lattice parameter of the 
regular structure) and indicates the coordinate of the (a;, 0, |) 
position of the CI atoms. The upper axis shows the change 
in the Cu-Cl-Cu angle relevant for the interaction Ji (ip, see 
the middle panel of Fig. HJ. 



tion Ji of about 30 K or 50 K depending on the position 
of the CI atoms. Other in-layer hoppings are relatively 
small, and the second largest hopping is the intcrlayer 
term (t±) that yields the AFM intcrlayer interaction of 
about 15 K. There is a pronounced difference between the 
inequivalent fourth-neighbor interactions. One of these 
interactions is comparable to J± (J^ FM ~ 20 K), while 
others are negligible. The interactions beyond the fourth 
neighbors are also negligible (the respective hoppings are 
below 10 meV implying J^ FM below 1 K). Thus, the TB 
fit suggests the following model for the spin system of 
(CuCl)LaNb207: AFM chains run along the b direction, 
and these chains are coupled by the AFM interactions J4 
and J± (see the right panel of Fig. [4]). The evolution of 
the leading hoppings parameters along the shift of the CI 
atom is shown in Fig. [7] 

The TB estimates are in qualitative agreement with 
the experimental data. The largest AFM coupling does 
not exceed 50 K and shows the same order of magnitude 
as the C urie- Weiss temperature^ and the saturation field 
(see Sec. Ill I J This strongly supports the structural dis- 
tortion and the proposed spin model in contrast to the 
frustrated square lattice scenario suggested by the reg- 
ular structure. The sizable long-range interaction J4 is 
in agreement with the INS dataP^ In Sec. [v] we will 
further compare our spin model with the experimental 
findings. Now, we proceed to the LSDA+C/ results that 
yield supplementary information on the leading exchange 
couplings in (CuCl)LaNb20y. 

For the LSDA+J7 calculations, we select all the in- 
layer interactions exceeding 1 K, i.e., Ji, J[, J2, J3, and 
J^J'^). To evaluate these couplings, one has to use a 
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TABLE II: Leading hopping parameters (U) of the tight-binding model and the resulting antiferromagnetic contributions to 
the exchange couplings (J* FM = 4t%/U e s) for two different positions of the CI atoms: x = 0.136 (experimental) and x = 0.16 
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TABLE III: LSDA+t/ estimates of the exchange couplings 
Ji (in K) in the distorted structure of (CuCl)LaNb207 (the 
positions of the CI atoms are fixed at x = 0.16). Ud (in eV) 
is the Coulomb repulsion parameter of LSDA+(7. 
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large, 96-atom supercell that makes full-potential calcu- 
lations quite time-consuming and, likely, not very accu- 
rate. Therefore, we employ an alternative approach and 
consider smaller cells. Three 48-atom supercells (Fig. [8| 
enable to evaluate a number of linear combinations of 
J\, J[, J%, J3, and (J4 + J4). Then the individual cou- 
plings are evaluated by solving a simple system of linear 
equations. We also use one 72-atom supercell (Fig. [8]) to 
distinguish between J4 and J4. The calculations for four 
different supercells are basically redundant, and some 
quantities are independently calculated within two dif- 
ferent supercells. The respective estimates match within 
2 K and indicate the internal consistency and complete- 
ness of our LSDA+[7 results. 

The LSDA+C/ estimates of individual exchange cou- 
plings are listed in Table |III| The calculations were done 
for x = 0.16 only. The results are in good agreement 
with the TB estimates. The leading AFM interaction is 
Ji, whereas the interactions J[ and J3 are weak consis- 
tent with the small hoppings t[ and t 3 , and J4 is dras- 
tically different from J' A . Ferromagnetic contributions to 
the exchange couplings are small. The values of J2 and 
J4 depend on the specific choice of the Ud parameter: at 
Ud = 7.5 eV, Jo is heavily overestimated with respect 
to J 2 AFM (Tablc"jnJ). On the other hand, the low U d of 
3.5 eV leads to a somewhat overestimated J4. Thus, the 
intermediate Ud value of 5.5 eV seems to be the optimal 
one, suggesting similar AFM interactions J2 and J4 that 
amount to approximately half of J\ . Basically, the choice 
of Ud — 5.5 eV also looks reasonable from the microscopic 
point of view. The Ud values of 6.5 — 8.0 eV are normally 
used for the LSDA+C/ calculations of Cu +2 -containing 
oxides with copper surrounded by oxygen atoms only, 




FIG. 8: (Color online) Different unit cells and supercells em- 
ployed in the present work. Small spheres show the Cu atoms. 
The dotted line indicates the 12-atom unit cell of the regular 
(tetragonal) structure of (CuCl)LaNb207. The dash-dotted 
line shows the 24-atom unit cell of the distorted structure. 
The dashed and solid lines denote the 48- and 72-atom super- 
cells employed in the LSDA+C/ calculations. 



while CI 3p orbitals provide better screening of the on- 
site repulsion thus reducing Ud^ 

The resulting exchange couplings can be rationalized 
as follows. In the distorted structure, the basic structural 
element is the CUO2CI2 plaquette. The strongest inter- 
action Ji runs between the corner-sharing plaquettes via 
the Cu-Cl-Cu superexchange pathway. The magnitude 
of this interaction depends on the respective hopping, 
while the hopping is controlled by the Cu-Cl-Cu angle 
((p, see Fig. [4]) and the position of the CI atom. In the 
regular structure, ip — 90°, the hopping t\ r is close to 
zero, and the nearest-neighbor interaction is weak FM 
(see Sec. III). The shift of the CI atoms leads to the in- 



crease in the ip angle and to the increase in the hopping 
t\ (see Fig. [7]). Then, the sign of J\ is changed and the 
interaction becomes AFM. At sufficiently large distor- 
tion, Ji is the leading AFM interaction. Yet, Ji in the 
distorted structure is one order of magnitude below J-2 r 
in the regular structure, since the ip angle for Ji is still 
close to 90° (e.g., <p — 111.5° at x = 0.16), while Ji r cor- 
responds to the linear Cu-Cl-Cu superexchange pathway 
(</j= 180°). 
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a crossover regime at weak distortions. In fact, the sce- 
nario of the small distortions is even more complex, be- 
cause the reduction in x leads to the enhanced fifth neigh- 
bor interaction (at x — 0.12, t& is about 30 meV). This 
scenario may be relevant for Br-containing compounds 
that show smaller distortions as compared to the Cl- 
containing counterparts. The respective discussion lies 
beyond the scope of the present work and should be a 
subject of a further study. 



V. DISCUSSION 



FIG. 9: Plots of the Wanner functions for the distorted struc- 
ture of (CuCl)LaNb2C>7 with x = 0.16. Left panel: the single 
CUO2CI2 plaquette. Right panel: the Cu-Cl-Cl-Cu superex- 
change pathway causing the large fourth neighbor interac- 
tion J4. 



In contrast to t\, the shift of the CI atoms does not 
change t[, because the respective Cu-Cl-Cu superex- 
change pathway includes one long, non-bonding Cu-Cl 
distance, implying the very weak overlap of the Cu and 
CI orbitals. The small hopping t\ suggests a possible 
FM interaction J[. The next-nearest-neighbor hopping 
is also small due to the lack of the proper superexchangc 
pathway. This interaction is reduced with the increase 
in the tp angle. Yet, some long pathways (about 9 A 
and 12 A for J4 and Jj_, respectively) are quite efficient. 
The origin of these long-range couplings can be under- 
stood by considering the Wannier functions centered on 
copper sites (see Fig. 9j. The single Wannier function 
is formed by the in-plaquette d orbital of copper with 
large contributions from tr-type orbitals of oxygen and 
chlorine (2p and 3p, respectively) . The oxygen p orbitals 
further overlap with niobium 4c? orbitals, hence the inter- 
layer coupling J± should involve the complex Cu-O-Nb- 
O-Nb-O-Cu superexchange pathway. In case of J4, the 
interaction is likely mediated by the direct Cl-Cl contact, 
as indicated by the overlap of the Wannier functions in 
the right panel of Fig. [9] This regime can be understood 
as the cr-overlap of the Wannier functions, because their 
lobes point to each other. The similar coupling regime 
has been reported for the Cu 2 Tc 2 5 X 2 (X = CI, Br) 35 
compounds that show sizable long-range interactions me- 
diated by the X-X contacts (the respective hoppings are 
about 80 meV). In contrast, the pathways J' 4 and J" do 
not lead to the effective overlap of the Wannier function 
and correspond to the very weak exchange. The same 
holds for another tellurium-containing copper oxychlo- 
ride Cu 4 Te 5 0i2Cl4 (Ref. d) where the Cl-Cl contact 
leads to the small hopping td < 10 meV. 

It is worth to note that the long-range interaction J4 
is basically independent on the precise CI position within 
the studied range of x = 0.12 - 0.20 (see Fig. [7]). On the 
other hand, the regular structure does not lead to any siz- 
able fourth neighbor interactions. Thus, there should be 



The spin gap behavior and the lack of the long-range 
magnetic ordering in (CuCl)LaNb 2 07 have remained a 
puzzle since the first report of Kageyama et aZP^ ap- 
peared. Further experimental studies^ indicated a struc- 
tural distortion, but the distorted str ucture remained un- 
known. Two models were suggested* 13 ^ 17 ! however, nei- 
ther of them was supported by total energy calculations 
and quantitative analysis of the exchange couplings. Our 
study presents a different structural model which is pro- 
posed and justified by careful full-potential band struc- 
ture calculations. We argue that this model is able to ex- 
plain the unusual magnetic behavior of (CuCl)LaNb207. 

Structural distortion is a key to understanding the 
properties of (CuCl)LaNb207. The distortion is realized 
via the cooperative shifts of the CI atoms. The shifts lead 
to the formation of short Cu-Cl bonds and the resulting 
CUO2CI2 plaquettes. The plaquettes are arranged per- 
pendicular to the [CuCl] layers and share corners to form 
chains (Fig. EJ. The origin of the distortion deals with 
the tendency to lift the orbital degeneracy of copper. The 
regular, tetragonal structure of (CuCl)LaNb 2 07 leads to 
two competing copper orbitals near the Fermi level, while 
the distortion stabilizes the single-orbital ground state 
with the half-filled in-plaquette orbital (Fig. [6]). Then, 
the resulting physics is similar to that of conventional 
cuprates with CUO4 square-planar structural units. The 
overall behavior resembles copper compounds with the 
cooperative Jahn- Teller distortion (e.g., KCUF3). How- 
ever, the magnitude of the distortion is different, and 
the driving force is likely different as well. In KCuF 3 , 
the shifts of the ligands are small (below 0.2 A), while in 
(CuCl)LaNb207 the CI atoms are shifted for about 0.6 A. 
To stabilize the distort ed str ucture of KCUF3, one has to 
include the correlations 1121201 wn ji e ^ ne distorted structure 
of (CuCl)LaNb 2 7 is stable even in LDA (Fig.[5|. 

The structural distortion has a pronounced effect on 
the magnetic interactions. The regular structure corre- 
sponds to the two-dimensional, frustrated square lattice 
model, while the distortion leads to a more complex spin 
model with numerous inequivalcnt exchange couplings. 
This model can be understood as one-dimensional in case 
only the leading coupling J\ is considered. Yet, the in- 
troduction of the interchain couplings leads to a three- 
dimensional spin model, since J^, and J± are of simi- 
lar magnitude. Band structure calculations yield numer- 
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ical estimates for individual exchange couplings, and it 
is instructive to compare these estimates with the exper- 
imental data. 

In the previous section, we have emphasized that the 
calculated J values have the same order of magnitude as 
the relevant experimental numbers, the Curie- Weiss tem- 
perature (9) and the saturation field ({j, H s ). To make a 
numerical comparison, we employ the spin model shown 
in the right panel of Fig. [4] The Curie- Weiss temper- 
ature is the coefficient at the second-order term of the 
high-temperature expansion for the magnetic suscepti- 
bility (x in powers of 1/T). Thus, 9 can be expressed as 
9 = 1/4 Y]j ZiJi, where Zi is the coordination number for 
the interaction Ji. The 9 value of 9.6 KP^is below the po- 
sition of the susceptibility maximum (T^ ax — 16 K) and 
the spin gap (A ~ 26 K)p3 Then the low Curie- Weiss 
temperature is likely caused by the presence of both FM 
and AFM interactions in (CuCl)LaNb207. Our calcula- 
tions do not show any sizable FM interactions and yield 9 
of 40—50 K, well above the experimental value of 9.6 K.H21 
Thus, the FM interactions are underestimated. Yet the 
AFM interactions seem to be overestimated, because the 
leading AFM interaction J\ ~ 50 K is above the satura- 
tion field of 30 T (i.e., 40 K)P However, the J values 
strongly depend on the position of the CI atoms, and one 
may expect a better agreement in case the experimental 
position is considered. Indeed, at x = 0.136 the leading 
coupling J i is reduced to 30 K (see Table |ll| , i.e., below 
the value of the saturation field. 

Basically, one has to find out the precise position of the 
CI atoms in order to give accurate estimates of the ex- 
change couplings in (CuCl)LaNb207. Nevertheless, our 
spin model derived for the optimized position (x = 0.16) 
is qualitatively valid and suggests the plausible mecha- 
nism of the magnetic frustration in this material. The 
right panel of Fig. [4] presents the spin lattice with the 
leading AFM interactions Ji, J2, and J4. The respective 
bonds form triangles; therefore, the three interactions can 
not be satisfied simultaneously. According to the second 



line of Table III J2 and J4 are of similar magnitude, 
hence a strong frustration of the interchain couplings is 
expected. Note that the FM interaction J[ also leads to 
a frustration via the four-site triangles formed by one J{, 
one J a, and two J\ bonds (i.e., one FM and three AFM 
bonds) . 

The strong frustration naturally explains the lack of 
the long-range ordering in (CuCl)LaNb 2 07 despite siz- 
able couplings along both a, b, and c axes being present. 
The explanation of the spin gap behavior is more com- 
plicated, because the spin model presented in Fig. |4|has 
not been studied theoretically. Whangbo and DaP^ re- 
ported the classical treatment of a simplified model with 
Ji and J4 only. They did not find the spin gap, and 
this result is not surprising, since the J\ — J4 model is 
non- frustrated. Normally, two-dimensional spin systems 
without the dimerization arc gaplcss. However, the mag- 
netic frustration can strongly modify the physics and l ead 
to a sizable spin gap as found in, e.g., CaV 4 9 PII3] We 



believe that a detailed theoretical treatment of the real- 
istic spin model with at least three interactions (Ji, J2, 
and J4) is necessary to elucidate the origin of the spin 
gap in (CuCl)LaNb207. Meanwhile, we should note the 
"unconventional" origin of the spin gap. The gap is not 
caused by dimerization; therefore, the gap does not sim- 
ply separate singlet and triplet states of a dimer. Rather, 
it separates a singlet ground state and a complex excited 
state of a two-dimensional frustrated spin system. The 
latter result may be relevant for the explanation of the 
drastic difference between the s pin ga p estimates from 
the high-field magnetization dat a 1 11 * 16 -* and the magnetic 
susceptibility or INS dataPl 

Thus, our calculations suggest a valid spin model that 
accounts for most of the experimental observations re- 
ported so far. The spin model is based on a structural 
model which is consistent with the available diffraction 
data and with most of the NMR results. According to 
Sec.|rV there is an additional distortion in the [LaNb207] 



block, which for now remains unresolved computation- 
ally. However, there is another structural aspect that 
may also be important. The calculations propose the en- 
ergetically favorable structural model shown in Fig. [4] 
This structure is thermodynamically stable as well, in 
case the entropy contribution is negligible (this is usually 
the case in ordered solids at sufficiently low temperatures, 
e.g., at room temperature). Nevertheless, the compound 
itself is metastable and readily decomposes at 450 °C. 
The preparation temperature is quite low (325 °C); there- 
fore, the real system does not necessarily reach the equi- 
librium state proposed by the calculations. Although we 
expect that the real structure resembles the thermody- 
namically stable one, the difference between the calcu- 
lated (equilibrium) and real (non-equilibrium) structures 
along with the uncertainty of the CI atom position are 
the most likely reasons for the discrepancies between the 
experimental and computational results. 

Unfortunately, there is no way to optimize the struc- 
ture of non-equilibrium systems. Therefore, our spin 
model and our estimates for the exchange couplings 
should be considered as approximations to the proper- 
ties of the real material. To get further insight, we 
suggest additional experimental studies. First of all, 
careful electron microscopy investigations are able to 
provide supplementary information on the structure of 
(CuCl)LaNb207. These studies should challenge the pro- 
posed structural model and investigate the local struc- 
ture, since a number of defects (e.g., twinning) are ex- 
pected. It is also desirable to use the experimental data 
to quantify the spin model proposed by our calculations. 
The Curie- Weiss temperature and the saturation field 
hold a certain quantitative information on the exchange 
couplings. Yet, the INS data should be more informa- 
tive. In Ref. Uni the INS data are analyzed within a 
simple dimer model. Further analysis within the realistic 
spin model (Fig. [1]) will be very helpful. 

Finally, we suggest that the scenario of the structural 
distortion may be generic for the whole family of lay- 
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ered perovskite-type copper compounds [(CuX)LaM 2 C>7, 
(CuX)M' 2 M 3 Oio, etc. with X = CI, Br; M = Nb, Ta; 
and M' = Ca, Sr]. All these materials present an in- 
teresting but presently poorly understood physics. For 
example, (CuBr)LaNb 2 Oy reveals columnar AFM or- 
dering at the unexpectedly high temperature of 32 K 
(despite the exchange couplings are also of the order 
of 30 - 40 K)P while some of the (CuBr)M' 2 M 3 Oio 
compound s sho w a magnetization plateau at 1/3 of the 
saturation.ESEfil The common feature of this compound 
family is the high-symmetry environment of copper and 
the tendency of the X atom to be shifted away from 
the high-symmetry position, as evidenced by high Debye- 
Waller factors in crystal structure refinements! 8 * 9 * 46 ^ We 
propose that all the compounds should in fact show a 
structural distortion of Jahn-Teller type. Then, shorter 
Cu-X bonds will lead to peculiar magnetic interactions, 
including long-range couplings. The first indication for 
the relevance of this scenario is found in the recent 
NMR study of (CuBr)LaNb 2 7 P It is shown that the 
bromine compound lacks the tetragonal symmetry sim- 
ilar to (CuCl)LaNb 2 07. This may suggest a distortion 
pattern similar to that in the chlorine compound (Fig.EJ. 
However, the different size of bromine will lead to a dif- 
ferent position of the X atoms, hence drastically changing 
the exchange couplings and the resulting physics. Fur- 
ther computational and experimental studies of the prob- 
lem are desirable. 

In conclusion, we have presented a microscopic model 
for the structural distortion and the resulting magnetic 
interactions in (CuCl)LaNb 2 Oy. The model is based 



on ab-initio electronic structure calculations, conforms 
to the experimental data (x-ray and neutron diffraction, 
NMR, and NQR), and explains the peculiar magnetic 
properties of the compound. We argue that the tetrag- 
onal crystal structure leads to a competition of copper 
orbitals, while the structural distortion lifts this com- 
petition and yields a single-orbital ground state. The 
pronounced distortion is stabilized by both lattice and 
correlation effects. The distortion gives rise to a drastic 
change in the magnetic interactions in (CuCl)LaNb 2 07. 
The lack of the tetragonal symmetry and the long-range 
interactions lead to a complex spin model with chains 
running along one of the directions in the ab plane and 
sizable interchain couplings. The interchain interactions 
in the ab plane are frustrated and may account for the 
spin-gap behavior. 
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